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Abstract 

The evolution of a class of inhomogeneous spherically symmetric 
universe models possessing a varying cosmological term and a material 
fluid, with an adiabatic index either constant or not, is studied. 
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One of the main puzzles concerning our current understanding of the 
physical world is the present small value of the effective cosmological con- 
stant | A |< 10~ 120 Mp,, or, which amounts to the same thing, the small 
value of the vacuum energy density p v = A/(8n G), as witnessed by cosmic 
observation. The status of the problem was reviewed by Weinberg [Qj. A 
current of thought holds the view that the cosmological term is not really 
constant but its value decreses as the universe expands. The rationale be- 
hind this is that the energy of the vacuum should spontaneously decay into 
massive and massless particles, hence reducing A to a value compatible with 
astronomical constraints, see e.g. §,@, §,§,§, |],§, 0, 0, 0, 



[y]. All these papers adhere to one or another "law" for the decay 



of p v . However, as shown by Pavon [14| these laws should be restricted by 
statistical physics considerations, and so several of them may be ruled out. 
All these studies were carried out on the assumption that the universe is 
homogeneous and isotropic, which is to some extent very reasonable, for it 
is believed the universe has been so at least since shortly after the begin- 
ning of the expansion. To our knowledge the only work departing from a 
Robertson- Walker background is the one by Beesham who simultaneously 
considers A and G evolving with time in a Bianchi type-I universe [|15|J. 

In this short article we consider a varying cosmological term in a class 
of inhomogeneous (but isotropic) universe. The reason to undertake such 
a study is that there is no reason a priori to believe that the universe was 
homogenoeous and isotropic right back to the beginning of the expansion, 
and on the other hand there is the fact that the homogeneous and isotropic 
models cannot account for the high degree of homogeneity and isotropy so 
frequently adscribed to the present state of our universe. Further motiva- 
tions to study inhomegeneous cosmologies can be found in [|1| as well as in 



[17]. 



Let us consider the isotropic but inhomogeneous spherically-symmetric 
spacetime described by the plane Lemaitre-Tolman-Bondi metric pi 



ds 2 = -dt 2 + Y 2 dr 2 + Y 2 (d9 2 + sin 2 9d0 2 ), (Y = Y(r,t)) (1) 

(the "prime" denotes partial derivative with respect to the radial coordinate 
r), the source of the metric being a material perfect fluid of energy density 
p(r, t) and pressure P(r, t) plus the quantum vacuum. The corresponding 
stress-energy tensor reads 

T ab = (P + P)u a U b +{P- ^q) 9ab , (2) 
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where u a = <5* is the fluid fourvelocity. The nontrivial Einstein's equations 
take the form 

P-A = -^ f (Y*Yy, (4) 
Y (Y\ 2 Y' YY' _ 

y + lyl ~ Y 7 ~ YY 7 ~ °' ^ 

where the upper dot means partial derivative with respect to time and we 
have set 8irG = 1. Introducing the change of variables 

Y = f 2/3 

the last three equations become into 

4 / f 
4 f 

F-A = - 5 f (7) 

and 

/"/ /'/ I). (8) 

respectively. We incorporate to this system the widely-used equation of 
state for the material fluid, namely 

P = ( 7 - l)p (9) 

where very ofently the adiabatic index 7 is considered constant, though on 
physical grounds it may also depend on time. The latter possibility is very 
natural since if the quantum vacuum decays into a mixture of massive and 
massless particles, 7 must vary with time because both species of particles 
redshift at different rates. Equation (Jsj) readily implies 

f-F(t)f = 0, (10) 

where F(t) does not depend on the radial coordinate. Note in passing that, 
in the particular case, 7 = 1 equations (Q) and (|i~0|) lead to A = A(i). For 
7 / 1 we will have in general that A = A(r, t), though A = A(t) is also 
possible. 
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Introducing the factorization 

f(r,t) = R(r)T(t) 

in (|n]) we get 

f-FT = 0. (11) 

The latter has two independent solutions, and therefore the general solution 
to (0) is 

f(r,t) = R 1 T 1 (t)+R 2 T 2 (t). 
However we will restrict ourselves to the simpler case 

f(r,t)=R(r)T(t), (12) 

which is nonetheless fairly general. In this way (||) becomes an identity and 
equations (||) and ([?]) reduce to 



and 



A (13) 



4 T 

P = --- + A, (14) 

respectively. Note that because of the above factorization any regular func- 
tion R(r) will satisfy (|l0|). To get an equation for T(t) we substitute the 
right hand sides of (|l3|) and in @; it follows 

TT + (7 — 1) T 2 — ^ A T 2 = . (15) 

To solve it we resort to the change of variables 

T = Z n (n = constant) , (16) 
which leads to n = I/7 and 

Z-^AZ = 0. (17) 

Depending on the expression of A, which in what follows will be assumed 
position independent (the latter automatically implies that both p and P 
will depend on t only), different cases arise. 



4 



1. For 7 and A constants one obtains 



Z x = Ci cosh t + ip^j , (18) 

Z 2 = C 2 sink f t + ip 2 \ , (19) 



and therefore 



Y 1 = R^{r) cos0^ t + n), (20) 

Y 2 = R y\r) CT" sinh^ t + tpX (21) 



Note that solution (|20| ) does not present initial singularity, but solution 
(21) has a singularity at to = — • However both sets of solutions 
have a final inflationary stage. 

2. For 7 = constant and 

A 2 

A(i) = ^| (A§ = constant) , (22) 

the corresponding differential equation can be integrated by using the 
ansatz Z oc t m with m = constant. The general solution reads 

Z(t) = Ci t m + + C 2 t m - (23) 

where m± = \ =t v^+^o _ inflationary solutions may occur for large 
enough A 2 ,. 

3. For 7 = constant and 

A = X 2 t n - 2 (n^O, 2), (24) 

equation ( |17| ) becomes 

Z- ^ 7 A§t n - 2 Z = (25) 

and the general solution can be expressed as a combination of Bessel 
functions (see Ref. 18) 
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+ C 2 t^J_ 1/n (^^t^y (26) 

The behavior at the asymptotic limits depends on n. For < n < 2 
one has the following, (i) When t —* one obtains 

Z ~ d t + C 2 . (27) 

One can choose C 2 = to have the initial singularity at i = 0. (ii) 
When t — > oo there follows 

Z ~ t2-T cost n/2 . (28) 

For n < one has the following: (i) when t — > one obtains Z ~ 
t2~4 cos {t n l 2 + if) . (ii) When t —* oo one obtains Z ~ t . 

4. For 7 = constant and 



A(t) = X z + ce" at (29) 

where Aq, c and a are constants (with c < for mathematical conve- 
nience), equation ( |l7| ) becomes 

Z - ^7 [A 2 , + ce~ Q *] Z = (30) 

and the general solution can be expressed as a combination of Bessel 
functions (see Ref. 18) 



/37 

with 



' l^ C ^> (31) 



6 



in order to fix the initial singularity at t = 0. 

The asymptotic behavior near the initial singularity, when t — > 0, is 
given by 

Z ~ t. (33) 

At the final stage, when t — > oo and A — > Aq, one obtains the following 
asymptotic behavior 

Y ps J R 2/3 (r)e^^* . (34) 

Besides, from ( ^l|) we recover the same result in the far future. 

For the particular case Ag = the general solution of ( |30|) is given by 



Z = Cl Jo e^* j + C 2 Y j , (35) 

where lo is the Weber function of the second kind and zero order. 
In the limit t — > oo and A — > the final behaviour of the solutions, 
obtained from ( |33|) are 

Y^R 2 / 3 (r)t 2 ^ . (36) 
The same result can be easily obtained from (|17D by setting A = 0. 



For 7 = 7(i) and A = A(t) we shall find the expressions of both 
quantities and analyse the behavior of the solutions at late time. To 
do this we introduce the new varible s(t) in the following way 

T = T ef^ dt , (37) 



inserting (|37 ) into ([lq) we get 



s--^s + -A 7 s = 0. (38) 
7 4 
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This equation can be identified with 

1 



s + s n s + 



2n+l 



(n^-2), 



(39) 



(n + 2) 2 

which is reduced to a linear differential equation by making the sub 

n + 2 



stitution 1 16] 



obtaining 



n c\ + j ' v n cit ' 



v = 0, v(t) = c 2 + c 3 t, 



(40) 



(41) 



where c\, c 2 and C3 are arbitrary integration constants. Equations (|3 
and (|39| ) are the same if we define 



(42) 



and 



3 1 

4 7 (n + 2) 2 ' 



2n 



(43) 



Without loss of generality we choose C2 = —to, where to is some initial 
time, and C3 = 1. So, the last system of equations can be easily solved 
to obtain 



A(t) 



2n 



4C 2 (t-t ) 
37on 2 (n + l) 2 



1 + 



(t-t r +i 

c 



and 



lit) = 70 



1 + 



(t - to) 



n+l 



c 



(44) 



(45) 



where 70 and C are arbitrary integration constants. On the other hand, 
inserting (^Tj) in (f40|), the general solution of the nonlinear equation 
( [39] ) is found to be 



s(t) 



(n + l)(n + 2) (t-to) r 



C+(t-i )n+l 



1/n 



(46) 



Now, taking into account that at late time, t 3> to, we must have 
7 — ► 70, the restriction n < — 1 readily follows as can be seen from 
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(|45|). In addition the cosmological term vanishes in the same limit. 
Now, using this approximation we evaluate T(t) in (|37|), finding 



T(t) « T 
and therefore 

Y^R 2 / 3 (r)T^° 



(n+l)(n + 2) (t _ to) l 1 ^ 



■n 



(n+l)(n + 2) (t _ t I 2 /** 



(47) 



(48) 



It is worthy of note that, for t S> io we have both 7 — » 70 and A — > 0. 
So, using these limits in ( |2T|) we recover the solution given by (48). 



To investigate the singular structure of the plane Lemaitre-Tolman- 
Bondi metric (|l]), we calculate the curvature scalar by resorting to change 
of variables Y = f 2 ^ 3 used above 



n 



r\ f 1 ^ f f 

7 3/7 
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(49) 



and evaluate it at the points where the coefficients of the metric Y' 2 and/or 
Y 2 vanish. To do this we insert the Einstein equation (§) along with (|l~2|), 
© and (ID in @, obtaining 



ft 



4 . (2 
77A + 



7) Z 2 



z 2 



(50) 



Then we replace in the scalar curvature the expansion of A and the corre- 
sponding solutions near the point where they vanish. All the solutions we 
have found for 7 = constant, except (pop, have a singularity at t = 0, i.e. 
the big-bang singularity. 

In summary we have found the coefficients of the Lemaitre-Tolman-Bondi 
metric assuming that the early universe possessed a time varying cosmolog- 
ical term, and that the adiabatic index of the material fluid were either 
constant or not. (i) All the solutions we have derived contain an arbitrary 
function of the radial coordinate, (ii) For 7 = constant all the solutions, 
except (pp|) have a singularity at t = 0, i.e. the big-bang singularity, (iii) 
Constant as well as varying cosmological terms give rise asymptotically to 
exponential inflation -see (pc|), ( ^T|) and (j34|). (iv) For a varying cosmolog- 
ical term there exist solutions which behave as though the universe were 
asymptotically matter dominated at late times when 7 = 1 -see (|3f 
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None of the solutions found has a spatially-homogeneous limit for t — ► 
oo. This is so because no homogeneization mechanism, such as anisotropic 
pressures [pC|] , was assumed. Such a more general study will be undertaken 
soon. 

This work was partially supported by the Spanish Ministry of Educatio- 
nunder Grant PB94-0718, and the University of Buenos Aires under Grant 
EX-260. 
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